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AND SEMICONDUCTOR DEVICE MANUFACTURING METHOD EACH 
EMPLOYING BOUNDARY CONDITIONS 



5 CROSS REFERENCE TO RELATED APPLICATIONS 

This application is based upon and claims benefit of priority from the prior 
Japanese Patent Application No. P2000-402778 filed on December 28, 2000, the entire 
contents of which are incorporated by reference herein. 

10 BACKGROUND OF THE INVENTION 

1. Field of the Invention 

The present invention relates to a simulation method, a simulation program, 
and a semiconductor device manufacturing method. In particular, the present 
invention relates to a technique of simulating semiconductor device processes 
15 employing to boundary conditions. 

2. Description of the Related Art 

A simulation of the manufacturing processes or electric characteristics of a 
semiconductor device is usually restricted by various conditions such as a computation 

20 time and a computer memory capacity. To comply with the restrictions, the simulation 
is carried out by setting a finite calculation area. Especially, a simulation involving 
three-dimensional calculations needs a larger memory capacity and a longer 
computation elapse time than a simulation involving two-dimensional calculations, and 
therefore, is usually allocated with a narrow calculation area. Setting a finite 

25 calculation area for a simulation is achieved by setting boundary conditions on the 
boundary of the finite calculation area. 
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The boundary conditions influence a simulation result. This influence of the 
boundary conditions is often overlooked by a design engineer, and the design engineer 
is frequently unaware of a deviation in a simulation result from a true result to be 
derived from the simulation on the calculation area. When relatively simple boundary 
5 conditions such as mirror, fixed, or periodic boundary conditions are improper for a 
given simulation, it is difficult to choose proper boundary conditions for the simulation. 
If improper boundary conditions are set for the simulation, an unintended simulation 
result will be obtained. If the unintended simulation result is passed unnoticed to 
design a semiconductor device, it leads to a failure of the semiconductor device. 

10 

SUMMARY OF THE INVENTION 
An aspect of the present invention provides a method of carrying out a 
simulation with simulation data, including, determining whether or not the simulation 
data includes boundary conditions set for a boundary of a calculation area set for the 

15 simulation, computing the influence of the boundary conditions on the inside of the 
calculation area if the simulation data includes the boundary conditions, displaying the 
influence of the boundary conditions on the inside of the calculation area, prompting to 
enter an instruction whether or not the boundary conditions are changed, and if an 
instruction to make no change in the boundary conditions is entered, carrying out the 

20 simulation with the simulation data. 

Another aspect of the present invention provides a semiconductor device 
manufacturing method, comprising, designing a semiconductor device, outputting 
design data of the semiconductor device, simulating the design data of the 
semiconductor device employing a simulation data in connection with the 

25 semiconductor device, the simulating including, determining whether or not the 
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simulation data includes boundary conditions set for a boundary of a calculation area set 
for the simulation, computing the influence of the boundary conditions on the inside of 
the calculation area if the simulation data includes the boundary conditions, displaying 
the influence of the boundary conditions on the inside of the calculation area, prompting 
5 to enter an instruction whether or not the boundary conditions are changed, and if an 
instruction to make no change in the boundary conditions is entered, carrying out the 
simulation with the simulation data, and fabricating the semiconductor device according 
to the design data. 



10 BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a block diagram showing a simulation system according to an 
embodiment of the present invention; 

Fig. 2 is a flowchart showing a simulation method according to an embodiment 
of the present invention; 
15 Fig. 3A explains boundary conditions; 

Fig. 3B shows an example of a window provided by the simulation system of 
Fig. 1 achieving the method of Fig. 2; 

Fig. 4 shows another example of a window provided by the simulation system 
of Fig. 1 achieving the method of Fig. 2; 
20 Fig- 5 shows another example of a window provided by the simulation system 

of Fig. 1 achieving the method of Fig. 2; 

Figs. 6A, 6B, and 6C show the behavior of ions under various boundary 
conditions; 

Fig. 7 explains a simulation method according to an embodiment of the present 
25 invention; 
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Fig. 8 is a flowchart showing a semiconductor device manufacturing method 
according to an embodiment of the present invention; and 

Fig. 9 is a perspective view showing a simulation system according to an 
embodiment of the present invention. 

5 

DETAILED DESCRIPTION OF EMBODIMENTS 
Various embodiments of the present invention will be described with reference 
to the accompanying drawings. It is to be noted that the same or similar reference 
numerals are applied to the same or similar parts and elements throughout the drawings, 
10 and the description of the same or similar parts and elements will be omitted or 
simplified. 

The present invention is applicable to, for example, process simulations and 
device simulations carried out to design a semiconductor device. A term "simulation" 
includes process simulation and device simulation in this specification. Also an 
15 apparatus used to carry out the process simulations and device simulations is called a 
simulation system in this specification. 

(Simulation system) 

Figure 1 is a block diagram showing a simulation system 100 according to an 
embodiment of the present invention. The simulation system 100 has a CPU 111, a 

20 RAM 112, a memory 113, an input unit 114, an output unit 115, a boundary condition 
database 116, and a bus 117. Through the bus 117, the CPU 111 reads programs from 
the memory 113, to carry out various control operations according to instructions 
written in the programs. The RAM 112 secures a work area to temporarily store 
programs and data related to processes executed by the CPU 111. The memory 113 

25 stores, for example, a simulation program 113a and a boundary condition display 
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program 113b. The memory 113 may include magnetic, optical, and semiconductor 
storage media that are readable by the simulation system 100. Programs and data to be 
stored in such storage media may partly or wholly be transferred to the simulation 
system 100 from a remote location through a transmission medium such as an electronic 
5 network. 

Figure 9 is a perspective view showing an example of a computer system 170 
that materializes the simulation system 100. The system 170 has a main body 180. 
The main body 180 is connected to a display 171, a keyboard 176, and a mouse 181. 
The main body 180 incorporates an FD (floppy disk) drive 172 to read an FD 173 and a 

10 CD (compact disk) drive 174 to read a CD 175. The drive 174 may be a DVD (digital 
video disk) drive to read a DVD 175. The system 170 may be connected to an external 
drive 179 to read an external memory 177 or a tape 178 such as a DAT. The FD 173, 
CD 175, external memory 177, and tape 178 are computer-readable storage media 
storing simulation programs, etc., according to the present invention. The programs 

15 are installed in and executed by the system 170, to carry out simulations according to 
the present invention. 

Returning to Fig. 1, the input unit 114 is used to enter various data related to 
simulations into the simulation system 100. The received data are transferred to the 
CPU 111, RAM 112, memory 113, etc. The input unit 114 includes, for example, a 

20 ten-key pad, a keyboard, and a mouse. The output unit 115 is used to output data 
related to simulations. The output unit 115 includes, for example, a display and a 
printer. The database 116 stores various boundary conditions that are retrievable and 
editable. The boundary conditions stored in the database 116 are, for example, mirror 
boundary conditions, fixed boundary conditions, periodic boundary conditions, 

25 computation accuracies, computation times, and reference documents. A user may 



retrieve necessary data from the database 116 and set proper boundary conditions for a 
simulation to carry out. 

(Simulation method) 

Figure 2 is a flowchart showing a simulation method according to an 
5 embodiment of the present invention. This method is achieved by the simulation 
system of Fig. 1. In step S201, the simulation system receives simulation data using 
the input unit 114 (Fig. 1), to start a simulation. The CPU 111 reads the simulation 
program 113a and boundary condition display program 113b from the memory 113 and 
stores them in the RAM 112. According to the programs, the CPU 111 executes 

10 processes explained below. Although the programs 113a and 113b are independent of 
each in this embodiment, they may be integrated into one by, for example, including the 
codes of the program 113b in the program 113a. The simulation data received in step 
S201 relates to a semiconductor device to simulate and may include boundary 
conditions. Step S202 checks to see if the received simulation data includes boundary 

15 conditions. If boundary conditions are included, step S203 is carried out. 

If no boundary conditions are included in the simulation data, the CPU 111 
outputs, in step S210, a message indicating there are no boundary conditions. This 
message may be accompanied by recommended boundary conditions. Step S211 asks 
the user whether or not boundary conditions must be set. If the user determines that 

20 boundary conditions must be set, the user edits the simulation data and sets boundary 
conditions in step S212. Thereafter, step S203 is carried out. If the user determines 
that no boundary conditions are needed, step S209 carries out a simulation according to 
the received simulation data. 

According to the boundary conditions found in step S202 or set in step S212, 

25 step S203 executes the boundary condition display program 113b to generate real 



images in a calculation area set by the user and virtual images outside the calculation 
area and provides the generated images to the output unit 115. 

The operations carried out in step S203 will be explained in detail. To 
compute, for example, a wire-to-wire capacitance value, two-dimensional calculations 
5 in a wire crossing direction are carried out if a wire length is sufficiently long. Figure 
3 A explains boundary conditions set for a boundary that defines a calculation area 130 
to compute a wire-to-wire capacitance value. In the calculation area 130, 
cross-sections of two wires A and B are shown. Each wire has a side length L, and the 
wires A and B are spaced apart from each other by a distance S. The outside of the 

10 calculation area 130 is a non-calculation area. The calculation area 130 extends from 
the wire A by a sufficiently long distance D. To calculate the capacitance value, a 
small potential dV is applied to the wire B and a charge quantity Q generated in the wire 
A is measured. Then, the capacitance value C is calculated as C = Q / dV. Data 
related to the wires A and B, data related to the calculation area 130 such as a dielectric 

15 constant, and the boundary conditions on the boundary of the calculation area 130 are 
included in the simulation data received in step S201. 

The boundary conditions will be explained. The boundary conditions are set 
on a boundary and are used for a simulation. The boundary conditions are represented 
with, for example, values and expressions. The boundary conditions may be fixed, 

20 mirror, periodic, transmission, and infinite boundary conditions. The fixed boundary 
conditions set fixed values for a part or the whole of a boundary. The mirror boundary 
conditions make a boundary to reflect, for example, ions that reach the boundary. The 
periodic boundary conditions make, for example, ions moved out of a boundary to 
reenter from another boundary. The transmission boundary conditions allow, for 

25 example, ions to move out of and into a boundary. The transmission boundary 
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conditions may be classified into perfect transmission boundary conditions and partial 
transmission boundary conditions. The infinite boundary conditions set an infinite 
distance outside a boundary when considering the behavior of, for example, ions. The 
examples shown in Figs. 3A and 3B employ the mirror boundary conditions. 
5 In Fig. 3A, the wire-to-wire capacitance value may correctly be calculated if 

the distance D between the boundary and the wire A is sufficiently long. If the 
distance D is short, mirror images formed on the outside of the calculation area 130 due 
to the mirror boundary conditions affect the wires inside the calculation area 130, to 
hinder a correct computation of the capacitance value. This is because, when the small 

10 potential dV to calculate the capacitance value is applied to the wire B, the wire A is 
affected by all the wire B including a real image and virtual images. This problem is 
frequently overlooked by the user if the user is unfamiliar with simulations. 

To cope with this problem, step S203 displays the virtual images that appear 
outside the calculation area 130 due to the boundary conditions, as shown in Fig. 3B. 

15 The calculation area 130 of Fig. 3B corresponds to that of Fig. 3A. The calculation 
area 130 includes the real images of the wires A and B. The non-calculation area 
outside the calculation area 130 includes the virtual images that are formed according to 
the mirror boundary conditions. The mirror boundary conditions form a virtual wire A 
on the left side of the calculation area 130. On the left side of this virtual wire A, a 

20 virtual wire B is formed due to the real wire B. With these images, the user may easily 
understand whether or not the capacitance value to be computed is an intended one. 

In Fig. 3B, the embodiment displays the mirror images outside the calculation 
area 130, so that the user may easily grasp the influence of the mirror images. To 
attract more attention of the user, the wire to which a small potential is applied and the 

25 mirror images of this wire may conspicuously be colored or blinked. The simulation 
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system may extract useful information such as a computation accuracy and a 
computation speed related to the simulation from the database 116 according to the 
boundary conditions set by the user and display the extracted information in an advice 
window 131a shown in Fig. 3B. 
5 Figure 4 shows another example of a window provided by the simulation 

system 100. The user of the simulation system 100 is allowed to change the 
calculation area 130, and the simulation system 100 generates mirror images based on 
the changed calculation area 130, as shown in Fig. 4. To change the calculation area 
130, the user manipulates a pointer 131b with, for example, the mouse. When the user 

10 clicks an OK button 131c, the simulation system 100 rearranges the mirror images of 
the wires A and B. Instead of pushing the OK button 131c, the simulation system 100 
may rearrange virtual images in real time in response to the movement of the calculation 
area 130 manipulated by the user. 

After step S203 displays images outside the calculation area 130, the user 

15 refers to the displayed images and information and determines, in step S204, whether or 
not the boundary conditions must be changed. If no change is needed, the user clicks, 
for example, a simulation execution button with the mouse, to instruct the simulation 
system 100 to start a simulation. If any change is needed, the user clicks, for example, 
a boundary condition change button with the mouse. Then, the simulation system 100 

20 executes the next step according to the boundary condition display program 113b. 

In step S205, the simulation system 100 refers to the database 116 and 
generates virtual images outside the calculation area 130 according to boundary 
conditions other than those set by the user, and displays the generated images condition 
by condition as shown in Fig. 5. The images generated in step S205 are prepared by 

25 temporarily replacing codes expressing the user-entered boundary conditions with codes 



expressing other boundary conditions and by executing computations with the replaced 
codes. The virtual images in Fig. 5 may conspicuously be colored or blinked like the 
example of Fig. 3B. Advise windows 131a in Fig. 5 are used to display information 
related to the boundary conditions. 

Boundary conditions used to simulate an ion implantation process in 
semiconductor device manufacturing will be explained. The ion implantation process 
implants impurities such as ions into, for example, a silicon substrate, an oxide film, or 
a polysilicon layer. A simulation of the ion implantation process calculates a 
distribution of ions in a calculation area of a target material made of, for example, 
silicon or polysilicon according to the acceleration energy, dose, and implantation angle 
of ions. In practice, ions are emitted from an ion emitting apparatus toward a target 
such as a wafer. In a simulation, an area in which an ion distribution is calculated is 
set as a calculation area, and the quantity of ions introduced into a material in the 
calculation area is calculated according to the paths of ions entering the calculation area. 
An ion implantation simulation on a two-dimensional structure will be explained with 
reference to Figs. 6A to 7. 

Figure 6A shows a rectangular calculation area and ions entering the area from 
above the area. The top side of the calculation area serves as a start point of each ion 
path, and a distribution of ions will correctly be simulated with the left and right sides of 
the calculation area being provided with any of the mirror, periodic, and transmission 
boundary conditions. 

Figure 6B shows a rectangular calculation area and ions obliquely entering the 
area. The top side of the calculation area serves as a start point of each ion path. If 
the left and right sides of the calculation area are provided with the mirror boundary 
conditions, ions that enter the right side of the calculation area will not be counted in 
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computation. On the other hand, ions that hit the left side of the calculation area are 
reflected back into the calculation area, to form a part whose ion concentration is higher 
than a real concentration. 

Figure 6C shows a rectangular calculation area and ions obliquely entering the 
5 area. If the left and right sides of the calculation area are provided with the perfect 
transmission boundary conditions, no ions on the right side of the calculation area will 
be considered in computation. If the periodic boundary conditions are applied to Fig. 
6C, ions that exit from the left side of the calculation area will be reintroduced into the 
right side thereof at the same coordinate heights, to provide a correct computation result. 
10 The periodic boundary conditions, however, are applicable only to a periodic structure. 

To cope with these problems, the simulation system 100 expands a calculation 
area as shown in Fig. 7, if the calculation area is improper for the relatively simple 
boundary conditions such as the mirror, fixed, and periodic boundary conditions. For 
the left and right parts of the expanded calculation area, the simulation system 100 sets 
15 boundary conditions that involve simple calculations. Namely, the simulation system 
100 generates a new calculation area based on an original calculation area, to correctly 
calculate an external influence on the original calculation area, and displays the new 
calculation area for the user. In Fig. 7, the original calculation area set for, for example, 
an ion implantation simulation is expanded by a length L on each side. The length L of 
20 each expanded area is expressed as follows: 

L=l*tan? 

where "1" (a lowercase el) is the height of the original calculation area and "?" is an ion 
incident angle. If an ion penetrating depth is known, the length L is expressed as 
follows: 

25 L=l*tan?*k + m 
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where k is the ratio of the ion entering depth to the height "1" and "m" is a margin. 
These expressions are useful to properly determine an expansion area, to minimize an 
increase in a computation time due to the expanded calculation area. 

These techniques according to the embodiment of the present invention draw 
the attention of the user to the boundary conditions. The user may easily understand 
whether or not an intended simulation result is obtainable. The simulation system 100 
may display not only the expanded calculation area but also the influence of the 
expanded area on the original calculation area. Although this embodiment relates to 
two-dimensional simulations, the present invention is also applicable to 
three-dimensional simulations. 

In this way, step S205 displays virtual images outside the calculation area 130 
according to boundary conditions other than those specified by the user. In step S206, 
the user determines whether or not required boundary conditions are included in the 
displayed boundary conditions. If there are, the user employs, for example, the mouse 
to specify the window in which the required boundary conditions are displayed. 
According to the specified boundary conditions, the simulation system 100 updates the 
simulation data in step S208. According to the updated simulation data, the simulation 
system 100 carries out a simulation in step S209 and provides a simulation result to the 
output unit 115. The user studies the simulation result and carries out various 
processes including a semiconductor device designing process accordingly. If there 
are no required boundary conditions in step S206, the user manually edits the simulation 
data in step S207 and returns to step S201. 

This embodiment asks the user to specify one of the windows displaying 
various boundary conditions. Instead, the user may specify the four sides of a 
calculation area and manually set boundary conditions on the four sides, respectively. 

12 



This enables the user to set required boundary conditions by him or herself. The user 
may partly set boundary conditions on the sides of a calculation area. 

Step S205 may employ a Monte Carlo simulation technique for simulating the 
behavior of impurities such as ions. This will be explained. The Monte Carlo 
simulation handles implanted ions as particles and calculates the hitting, scattering, and 
stopping of the ions with respect to atoms in a target material into which the ions are 
implanted. In such calculations, ions that enter a calculation area scatter, change 
directions, and sometimes move out of the calculation area. If calculations are 
continued outside the calculation area, some ions exited from the calculation area will 
never return into the calculation area and some will be redirected into the calculation 
area. The Monte Carlo simulation is capable of correctly perform these calculations 
even if the target material has an intricate structure. A computation time of the Monte 
Carlo simulation extends in proportion to the number of particles in the calculation area. 
If the calculation area is unnecessarily expanded, a very long computation time is 
needed because a large number of particles will be included in the calculation area. 

To cope with this problem, the present invention preliminarily introduces a 
small number of ions into an expanded calculation area and examines positions where 
the ions stop. Based on the ion stopped positions, the present invention determines a 
final expanded calculation area and achieves a Monte Carlo simulation. At this time, a 
result of the computation made on the small number of preliminarily implanted ions is 
included in the Monte Carlo simulation, to save a calculation time. For example, about 
one percent of particles (ions) used for the Monte Carlo simulation may preliminarily be 
introduced to determine the final expanded calculation area. Only by increasing a total 
computation time by one percent, this technique provides a correct simulation result. 

In this way, the embodiment studies the influence of an expanded calculation 
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area on an original calculation area when simulating an ion implantation process by the 
Monte Carlo simulation technique. This eliminates simulation errors that may involve 
repetitions of calculations, and correctly simulates an impurity distribution. 

The embodiment mentioned above expands an original calculation area, 
5 displays an expanded calculation area, and asks the user to determine a final calculation 
area. According to an embodiment of the present invention, the simulation system 100 
may automatically expand an original calculation area, count the number of particles 
that exit from and return to the expanded calculation area, carry out statistical 
calculations, and according to a result of the statistical calculations, check to see if the 

10 influence of the expanded calculation area on the original calculation area is below a 
predetermined level. If the influence is below the predetermined level, the simulation 
system 100 starts a final simulation. This is a technique to automatically and 
efficiently perform a simulation. 

The embodiments mentioned above locally input and output data. According 

15 to an embodiment of the present invention, the data input and output may be achieved 
between remote locations through an electronic network. For example, a LAN may be 
employed to input and output simulation data, simulation results, and boundary 
conditions. 

As explained above, the simulation system 100 computes and displays a finite 
20 calculation area used for a simulation and the influence of boundary conditions on a 
simulation result. With the simulation system 100, even a user unfamiliar with 
simulations can determine boundary conditions suitable for simulating physical 
quantities related to a semiconductor device to develop. 

The simulation system 100 allows the user to properly change simulation data 
25 while examining the influence of boundary conditions on a calculation area and the 
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effectiveness of the calculation area. With the simulation system 100, the user can 
efficiently execute simulations. When applied to process and device simulations for 
designing a semiconductor device, the simulation system 100 greatly reduces labor and 
time needed for the simulations. With the simulation system 100, the user can 
examine the influence of boundary conditions on a simulation, and according to a result 
of the examination, can correct simulation data to execute the simulation. With the 
simulation system 100, the user can correctly evaluate boundary conditions set for a 
simulation and can set a proper calculation area that provides a correct simulation result. 

The simulation program according to the embodiment is applicable to process 
and device simulations in designing a semiconductor device. With the program, the 
user can properly evaluate boundary conditions set for the simulations and correctly 
design the semiconductor device. The simulation program greatly reduces the labor 
and time of semiconductor device designing. The simulation program enables even a 
user unfamiliar with simulations to easily and correctly carry out simulations with 
proper boundary conditions. 

(Semiconductor device manufacturing method) 

Figure 8 is a flowchart showing a semiconductor device manufacturing method 
according to an embodiment of the present invention. Step S301 is a design process to 
provide design data 302 related to the semiconductor device. According to the design 
data 302, step S303 carries out a simulation to evaluate the electric characteristics and 
manufacturing processes of the semiconductor device and determine physical quantities 
concerning the semiconductor device. The simulation is carried out to improve the 
characteristics of the semiconductor device. The simulation employs the simulation 
system 100 mentioned above, to compute the influence of boundary conditions 
contained in simulation data on a calculation area set for the simulation. According to 
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the computed influence of boundary conditions, final boundary conditions are 
determined. According to the final boundary conditions, a simulation is carried out to 
determine physical quantities concerning the semiconductor device. According to the 
physical quantities, step S304 manufactures the semiconductor device. In this way, the 
5 semiconductor device manufacturing method according to the embodiment involves the 
simulation process that employs the simulation system 100 of the present invention to 
easily and correctly carry out a simulation with boundary conditions. The 
semiconductor device thus formed has improved characteristics and precision. 

The present invention may be embodied in other specific forms without 
10 departing from the spirit or essential characteristics thereof. The embodiments are 
therefore to be considered in all respects as illustrative and not restrictive, the scope of 
the present invention being indicated by the appended claims rather than by the 
foregoing description, and all changes which come within the meaning and range of 
equivalency of the claims are therefore intended to be embraced therein. 
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